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Abstract 

We calculate the ss condensate in the D - and B - mesons using the unquenched 
hopping parameter expansion at infinite coupling for the Wilson lattice action. We 
discuss the phenomenological relevance of our result. 



1 



1 Introduction 



Sea quark effects have over the last few year become increasingly more important in ex- 
plaining a variety of experimental results. Examples are the proton spin as measured by 
the EMC |||, violations of the Gottfried sum rule obtained by the NMC [0] and confirmed 
by the NA51 Collaboration H as well as abundant c^-meson production in pp annihilation 
Q at LEAR. The experimental indications for sea quark effects are supported by recent 
lattice QCD calculations [||, || of the pion-nucleon sigma term where it was found that the 
sea contribution is twice the valence one yielding a value for the sigma term in agreement 
with the result from chiral perturbation theory Q. It is therefore reasonable to expect 
that sea quark contributions are important in a number of other observables. In particular 
they may have implications for hadronic decays of the D- and B- mesons. In ref. || it 
was pointed out that a significant amount of ss in the D- and B- mesons would enhance 
decays of these mesons to strange final states, thus providing an explanation for the large 
branching ratios || for D° — » cpK and D + — > (pK + . Similar implications would follow for 
the decays B — > <j) + X . 

In this work we estimate the ss content of these mesons using the unquenched hopping 
parameter expansion of the Wilson lattice action at infinite coupling. Infinite coupling cor- 
responds to quark propagation in a background gluon field. Because of this simplification 
it is feasible to perform a convergent hopping parameter expansion analysis of some expec- 
tation values. In the quenched approximation one can define propagators in the random 
walk representation and then sum to infinite order in the hopping parameter |U], p!2[ . In 
spite of the fact that the analysis is done at infinite gauge coupling one still obtains the 
meson masses to a few percent accuracy [11]. In the unquenched case one can no longer 
define propagators as random walks and one has to truncate the expansion to some order 
and estimate the systematic errors due to the truncation |l3| . In our work we truncate the 
unquenched expansion at the 8th order. This still leads to meson masses within 10% of 
their experimental values. 

The quantity of interest here is the ratio 

< M\uu\M > 

< M\ss\M > ^ 

where M stands for a D- or a B- meson. It is shown in the next section that in the 
unquenched hopping parameter expansion this ratio can be obtained by differentiating the 
mesonic mass. We find that the ss contribution is suppressed as compared to the valence 
contribution by at least an order of magnitude. We proceed to explain how we reach this 
result. 



2 Calculation of mesonic expectation values 



Our starting point is the Wilson lattice action of QCD [14]: We consider an SU(3) color 
matrix U(b) in the fundamental representation defined on each oriented lattice bond b. 
Our convention is that 

U(-b) = C/t(6). (2) 
An oriented path lo on the lattice is a set of bonds 

w = bi U b 2 U . . . U b n (3) 
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such that the end-point of bi is the starting point of for 1 < i < n — 1. We can associate 
an SU(3) color matrix with u> by defining the path ordered product 

U(u) = U(b l )U(b 2 )...U{b n ) . (4) 

The spin matrices are defined in terms of 7-matrices by 



(5) 



Y(h) — 1 = 1 + 7^ if b in direction 
1 f v = 1 — 7^ if b in direction 

where the 7-matrices are hermitian 4x4 matrices, satisfying {7^,7^} = 25^. 
The Wilson action is then defined on a lattice A by 

S = -/3 Yl TrU(dp) - £ \k f ]T i; f (x)T(b)U(b)^ f (y)-Y^ f (x)^ f (x)\ (6) 

pCA* f=flavors [ b=(xy) x£A ) 

where p represents a plaquette and the kf are the hopping parameters associated with the 
different quark flavors. A* denotes the set of all bonds defined on A. The quark fields are 
represented by the anti-commuting variables ipf(x) and ^/(x) which transform under the 
3 and 3 representation of color. 

The mesonic expectation value of interest here is then obtained by evaluating 

(M^(x)\^ f ^ f3 \M^(y)) = 

= ±J [d##U]M^{x)^M^{y) W/s(*)) e- 5 ^) (7) 

where Z is the partition function and M^ 2 (x) = ip^ (x)Mipf 2 (x) represents a meson 
operator built up of an anti-quark of flavour f\ and a quark of flavour f 2 (M is some 
matrix in spin and color space). The three point functions that must be evaluated in order 
to obtain the desired ratio are shown in Fig. 1. They evolve a valence contribution when 
fs is equal to /1 or f 2 and a sea quark contribution. 

To avoid the calculation of three point functions one introduces a parameter A/ in the 
action 

S x = -P Yl TrU(dp)- Yl \ k f E MWWWfiy) - X f £ ^(x)^( 3 

pCA* f=flavors b=(xy) x&K 

(8) 

Since the meson propagator is given by 

(M /l/2 (x) t M /l/2 (y)) A = i- / [d^][d^][dU}M hh {x)^M hh {y)e- s ^ u ^ ) (9) 

A J 

the derivative of the meson propagator 



(M|^ / 3^ /3 |M) conn . = -JL(Mf^(x)lM^(y)) x 



(10) 

A/=1V/ 
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(a) 




(b) 

Figure 1: (a) Valence contribution and (b) sea quark contribution to the ratio given in (1) 



leads to the desired expectation value of the connected diagrams in ([?]) after setting this 
parameter to unity to recover the original action. 

The ratio between two different mesonic expectation values can be expressed as the 
ratio of the derivative with respect to the A's of the meson mass, assuming that for large 
space-time separation a meson propagator behaves proportionally to e~ mM ' x '^L 



C( 

(M|V7'V7'l M )c 



dm M (\)/d\ f 



dm,MW/dXf 



(11) 



A/=1V/ 



3 Hopping parameter expansion for the meson propagator 

Mesonic masses are obtained by studying the long distance behaviour of the propagators 
given in (^) . In this work the propagators are evaluated using the the hopping parameter 
expansion at the infinite coupling limit. This is done by breaking the action S into two 
parts, the unperturbed part Sq and the perturbation Si: 

5 = £^M E^/O^/^) f 

/ ixeA ) 

51 = -El*/ £ 4> } (x)T{b)U(b)^ f (x)\. (12) 

The plaquette term is not present because we take the infinite coupling limit i.e. f3 — 0. 
Expanding the exponential in eq. (S) in term of the perturbation and integrating out the 
gauge degrees of freedom U we obtain an expression for the the meson propagator 

(MM»(x)iM*Hv)h= £ C M A(u+Uu- n )k^W£ l "£~ k fL ( 13 ) 
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where 

A{uj + \Juj~ :1 ) = Tr (M t r(w + )Mr(w^)) Tr(T{j)) (-l) L -> (14) 

In ([[^) we have rescaled the hopping parameters kf to kf = kf/Xf. The term Cm is an 
overall constant amplitude which does not affect the meson mass, u + is a path on the lattice 
connecting x with y, u)~ is a path connecting y with x, uj = u) + U ui~ is the combination 
of them, 7 is a closed path representing the sea-quarks loop contributions and L 7 is the 
number of closed loops formed by the path 7. The mapping 0(o;, 7) denotes the group 
integral over the Haar measure of SU(3) associated with the graph defined by to and 7. 
Finally, T(oj) and U (to) represent the ordered products of T(b) and U(b) on the bonds b E uj, 



respectively. We do not give a proof of (J13, 14 ) since it is standard fll5| , [T 



In order to compute the meson masses we consider the static propagator 

G M (i) = £(M((0,0))tM((t,*))) (15) 

X 

(we omit the flavour indices for simplicity). The lowest order diagram representing a static 
propagator is a double fermion path uq in the time direction (Fig. 2a). The full propagator is 
given by a sequence of exited states connected by lowest order states shown schematically in 
Fig. 2b where each bubble can be viewed as a process contributing to the excitation. These 
excitations renormalize the mass of the static propagator. The unrenormalized meson mass 
is given by the lowest order diagram of the static propagator, namely 

Gi\t) = C M x exp(-mot) = C M x {Ak h k h f. (16) 

This tells us that the unrenormalized mass is niQ = —log(4:kf 1 kf 2 ). Let us now consider 
the effect on the mass produced by intermediate exited states. Fig. 2c shows an excitation, 
with the static propagator entering at the initial point w and departing at the final point 
z. For fixed w and z, we sum over all intermediate exited states to obtain the total weight 
for this event. We denote this weight by D M (w, z). The full propagator takes the form 

00 

(M(x)tM(y)} = E E G% I (x,w 1 )D M (w u z 1 )G% I (z 1 , W2 )D M (w 2 ,z 2 )x 

n=0 wi,Zi(i=l...n) 

x...xD M (w n ,z n )G^(z n ,y) (17) 

where the points Wi and Z{ are required to be time ordered and n represents the number of 
excitations. This last equation is represented pictorially in Fig. 2b. 



Following [15 1 one can factorize all static unrenormalized terms Go in eq ( |17|) and 



compute the excitations D(wi,Zi) separately. We define^ 

D M (t) = (c*Ht)r 1 Y, DM ^y) 



(18) 

yo-%o=t 



Notice that (G^(t))- 1 = ^e m ° ( = ^(4fc : 
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t 




where we normalize the excitation relative to the static unrenormalized propagator over the 
time interval t. Prom this definition and eq. we obtain 

G M (t) = e~ m ^ D M (t 1 -s 1 )...D M (t n -s n ) = 

Si,ti(i=l...n) 

e -m t e p M t ( 19 ) 



no 



where ti and Sj are the corresponding time coordinates of Z{ and W{. The exponential 



representation of (19) is up to an irrelevant boundary term exact and yields the renormalized 
mass m,M 

m M = m -p M - (20) 
The mass correction pu can be written in the form 

PM= E {G?V))- 1 D u {0,z)= £ (4A: 2 )- t 'D M (0,z) . (21) 

z=(t'>0,z) z=(t'>0,z) 



The last equation follows from ( |19| ) if we expand the exponential function containing pm- 
The excitation term D M (w, z) is defined starting from eq. (|i~3| ) in the following way: 



dM ^=a^cz £ ^#r^"'E^i (22) 

U oj:toi-^^^^;7e5(w) f s > 



./ .s C (I 



where the sum is over all closed paths u = cu + U u)~ from w to z and S'(w) C A* is the set 
of all loops 7 for with Q(uj, 7) 7^ 0. A' 1 -' 2 7) denotes the amplitude of the intermediate 
excitations, = 4 is a constant which has to be explicitly evaluated from (|l3|). 

To evaluate the correction to the unrenormalized mass mo due to the excitations from 
eq. (^TJ) we have to identify all closed paths to from the point to the point z = (t, z) 
and to sum over all points z. At (3 = we do not have plaquettes in the expansion 
because the plaquette expectation value vanishes, therefore a diagram consists of a closed 
line uj = lo + U u)~ composed by a set of connected bonds coming from the valence-quark 
contribution and a second path 7 coming from the sea-quark contribution. The group 
integral of a diagram is non- vanishing only if, for all b G co U 7, there exists a second link 
b 1 E WU7 such that b = —b' (b is equal b' in opposite direction). A method to easily evaluate 
the group integral is discussed in [|i~6[. 



4 Results and Conclusions 

We now evaluate the mass of the pseudoscalar and the vector mesons with fi and /2 quarks 
flavours up to order eight in the expansion, including valence and sea quarks. The sea quark 
contribution arises from quarks with the same flavour as the valence quarks, incorporated 
in the notation kfajkfy, as well as with flavours different from the valence quarks denoted 
by k sea . All terms in the expression of the form + kj^ + k™ e( ^j with n = 4 or n = 6 
correspond to sea quark contributions to the meson mass. Collecting all terms we obtain 
for the pseudoscalar meson mass 
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mp S ({k fl k f2 };k S ea) = (23) 

" ln(4 k h kf„ ) - 3 k h k h - 36 kl k% - 432 k% k% 



-19 § + kf e + kf ea ) - k h k h (kl + k% + kf ea 
-5328 $ - 8 q + kf 2 + kl a ) - 8 g (kl + kl + k 
-16 k h k h (k % h + kl + kL) ~ 32 k h k% (kl + k% + kf ea ) 
-32 § fc A (kl + kl + fc s 4 ea ) - 48 kl kl (kl + kl + ktea) 
+8*£ (kl+kl+kL)+8kl (kl+kl + k 4 s 
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sea 



and for the vector meson mass 



m VK ({k fl k f2 };k sea ) = (24) 

" M4 % k h ) ~ 2 fcfr % - 16 g g - 128 g § 
-3g(^+^+^ 4 ea ) 

-19 kl (kl + fc£ + ki a ) - k fl k h (q + kl + fc s 4 



"sen 



-1088 kl kl - 8 kl (kl + kl + kl a ) - 8 kl (kl + kl + *« 



-16 k fl k h (kl + $ + kl a ) - 24 % k% (kl + kl + kj 
-24 kl k f2 (kl + kl + fcl) - 32 kl kl (kl + % + fc s 4 ea ) 
+8 (kl + # + + 8 kl (kl + kl + 4 4 ea ) 



To fix the quark masses we use the experimental values of meson masses expressed in 
units of the p mass. For the u- , d- and s- quark masses we use the experimental pion 
and kaon masses (we consider the u- and the d- quarks to have the same mass). For this 
evaluation we neglect all sea quarks of flavour different from u, d and s. We thus obtain 
two equations 

,,, ., .in- — i- . i. ;. — 7. i 

(25) 



m PS ({k u = 


' kd\] k sea 


= k s ) 


_ m n 


m V K({k u = 


z k<i}', k sea 


= k s ) 


m p 


m PS ({k u , 


ks}'-, k sea - 


= 0) 


_ m K 


mvK({k u = 


: kd}', k S ea 


= k s ) 





~ 0.64 (26) 

from which we can extract the values of k u and k s . The resulting values are k u = kj = 
0.28(3) and k s = 0.25(2). The error given here is an estimate of the systematic error from 
the truncation of the expansion. It was obtained by comparing the truncated results to the 
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known bounds in the meson masses [O] ^| 

3k 2 2k 2 
mps < - log(4/c 2 ) - 1 _ 12fc2 myj( < - log(4/c 2 ) - (27) 

where k = max(/cj 1 , /c/ 2 , k sea ). 

Using these parameters one can check the expansion by predicting the K* mass obtained 
from the vector meson mass formula 

m VK ({k u = 0.28, k s =0.25}; k sea = 0) = 
m VK ({k u = k d = 0.28}; k sea = k s ) K ' { ' 

which compares favourably with the experimental value of mK*/mp — 1-16. Knowing the 
u- and s- quark masses one can estimate the ratio of the ss to uu condensates in the vaccum. 
To order eight it is easy to compute and we find < ss > / < uu >= (k s / k u ) 8 . Because this 
depends on the eighth power of the k— values it is very sensitive on the errors of k u and k s . 

In an analogous way one fixes the c-quark mass using the experimental value of the 
D- meson mass. The accuracy of the expansion can then be checked by comparing the 
predicted mass of the D* with the experimental value. In this case we find k c = 0.14(1) 
leading to a D* meson mass ratio mp>* /m p of 2.5(1) to be compared with the experimental 
value of 2.62. 

Having fixed the quark masses we obtain for the ratio 

(D\uu\D) dmp S (k u ; k c = 0.14; k sea = k s )/d\ u 



(D\ss\D) dmps{k u = 0.28; k c = 0.14; k s )/d\ s 



~ 50(10) (29) 

fc u =0.28;A; s =0.25 

The error given for the ratio is estimated by taking into account only the errors in the values 
of the A;'s which include part of the systematic error due to the truncation. The latter can 
not be estimated in the same way as the error in the A;'s because no upper bounds exist for 
the derivative of the meson mass. 

The same analysis can also be done in the case of the B- meson where one uses the 
experimental value of the B- meson mass to fix k^. We find kb = 0.011(2) leading to a B* 
meson mass ratio mg. /m p of 6.8(4) within the experimental value of 7.06 and for the ratio 

(Z?(/,,|Z?) 45(5) . (30) 



(B\ss\B) 

We note that, for the allowed parameter range (k u ,kj,k s ) G [0,0. 28] 3 where kf is 



the hopping parameter corresponding to the c- or the b- quark the ratios (|29 
bounded from below by 



remain 



<mHm) > 40 ^kf,k s ) e[0,0.28] 3 . (31) 

From the above values of the ratios we conclude that the non-valence contribution in 
the D- and B- mesons is suppressed compared to the valence contribution by at least an 



2 It was shown by the same author that a very good approximation to the meson mass is given by 
mps — — log [4kf 1 kf 2 /(l — 12/c/j kj 2 )] and mvK — — log [4fc/ x fc/ 2 /(l — 8fc / x fc/ 2 )]. The deviation of the 
truncated results from these values can also be used to estimate the systematic error due to the truncation. 
One obtains a similar value as the one quoted in the main text. 
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order of magnitude. This would comply with the OZI rule meaning that the OZI rule can 
not be evaded in the limit of infinite coupling. Therefore the enhancement observed in D- 
meson decays to strange final states seems not to be explained by a rich ss content in the 
D-meson. Whether finite (3 effects can change the above ratio by an order of magnitude 
can only be checked by doing a finite lattice simulation. 
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